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Abstract 

Hadron spectra and hadronic moments in the decay B — > X s £ + £~ are calculated taking into ac- 
count both the short-distance and long-distance contributions in the decay amplitude using a Fermi 
motion (FM) model to incorporate the B-meson wave-function effects. The measured branching 
ratios for the inclusive decays B — > X s + (J/tp, ip', ...) — > X s £ + £~ are used to fix the normalization 
of the long-distance contribution. The momentum distribution of the J/ip measured by the CLEO 
collaboration is fitted in the FM model which is then used to calculate the hadronic spectra from 
the resonant contribution also away from the J/tp- resonance. We also study the effect of various 
descriptions of the resonant and non-resonant cc contributions in B — > X s l + t~ existing in the 
literature on the hadron energy and invariant mass spectra, and in the Forward-Backward asym- 
metry. Selective cuts on the hadron and dilepton invariant masses can be used to reduce the BB 
background and resonant contribution and, as an example, we work out the hadron spectra with 
the experimental cuts used by the CLEO collaboration in searching for the decay B — > X s £ + £~. 
We show that data from the forthcoming B facilities could be used effectively to measure the short- 
distance contribution in B — > X s £ + £~, enabling precise determination of the FM model and heavy 
quark effective theory parameters Ai and A. 
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1 Introduction 



In two earlier papers HJ,^, we have worked out the short-distance (SD) contribution to the hadron 
energy and hadronic invariant mass spectra and the first two spectral moments in the inclusive decay 
B — > X s £ + £~ . The calculations reported in these papers were based on the leading order perturbative 
QCD corrections in a s and leading power corrections in using the heavy quark expansion tech- 

nique (HQET) (3|,|]]. In particular, we worked out the dependence of the moments and (Ejj), 

for n = 1,2, valid up to C(a s /m|, l/m|). It was argued that their measurements in forthcoming 
experiments could be combined with the improved measurements of the same in the semileptoic de- 
cays B — >• X£i>e to determine the HQET parameters Ax and A precisely. Since these parameters are 
endemic to most applications of HQET, their precise determination would reduce the present theoret- 
ical uncertainties improving the standard model calculations. In particular, the determinations of the 
CKM matrix elements V u b and would be considerably improved. The correlations resulting from 
(assumed) values of {Sh) an d {Sjj) were shown and compared with the constraints emerging from 
the analysis of the decay B — > X£vi reported in |j. The hadron spectra and spectral moments were 
also calculated in a phenomenological Fermi motion model (FM) 0. In particular, the remarkable 
similarity of the hadronic moments in the HQET and FM model approaches was quantified. Finally, 
the effects of the experimental cuts on the spectra used in the searches for the decay B — > X s £ + £~ by 
the CLEO collaboration |7j were studied. However, the effects of the long-distance contributions were 
not included in the hadron spectra or the spectral moments. 

The aim of this paper is to calculate the profile of hadron energy and invariant mass spectra 
by incorporating the effects of the long-distance (LD) contributions in the decay B — > X s £ + £~ . As 
opposed to the SD-contribution discussed in [Q,|2|, the LD-contributions are estimated phenomenolog- 
ically. To that end, the branching ratios for B — ► {J/ip,ip\ ...) + X s — > X s £ + £~ are described in the 
factorization approach (§|-|9| , with the data fixing the normalization and phase of the LD-contribution 
in B — ► X s £ + £~dX the J/ip,if)', ... resonances fllOfl . Assuming a Breit-Wigner form, one can extrapolate 
the dilepton mass spectra away from the J/ip,ip' , ... resonances. Using this and the SD-contribution, 
various distributions in B — » X s £ + £~h&ve been worked out in the literature [10-16]. We use the 
FM model to incorporate the wave-function effects. In the process of doing this, we also show that 
the FM model || provides an adequate description of the J/^-momentum distribution in the decay 
B — > J/ip + X s — » X s £ + £~ measured by the CLEO collaboration fl7|| . This confirms a similar and ear- 
lier study on this point |18|| . However, we have redone a fit of the CLEO-data on the J/^-momentum 
spectrum and prefer somewhat different FM model parameters than the ones presented in J18[| , mo- 
tivated by the analysis of the photon energy spectrum in the decay B — > X s + 7 |19| and theoretical 
consideration on the 6-quark mass |20t|. 
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The prescription of adding the SD and LD-contributions in the decay B — > X s £ + £~is not unique, 
which introduces a theoretical dispersion in the resulting spectra. Related to this is the inherent uncer- 
tainty concerning the extrapolation of the resonant part far away from the resonances 



study these uncertainties by working out a number of Ansatze used in the literature. The dispersion 
on the spectra emerging from various approaches can then be taken as a measure of theoretical sys- 
tematic errors from these sources. We analyze the dilepton invariant mass distribution, the Forward 
Backward (FB) asymmetry, hadron energy and hadronic invariant mass spectra and the hadron spec- 
tral moments in this context. The LD/SD-related uncertainties are found to be small for the hadron 
spectra and spectral moments. Some of these issues were also discussed in the context of exclusive 



decays B — ► [K, K*)£ + £~ in the second reference cited in 16]. 

Finally, we study the effect of experimental cuts on the branching ratios, hadron spectra and 
spectral moments in B — > X s £ + £~ . For that purpose we take the cuts used by the CLEO collaboration 
[f7j|, which involve dilepton and hadronic invariant masses. We present a comparative analysis of these 
cuts on the inclusive hadron energy and invariant hadronic mass distributions with and without the 
cc-resonant contributions. This shows that the cuts employed in |7| are effective in removing the 
resonant part and given data one could study the more interesting SD-contribution in B — > X s £ + £~ . 

This paper is organized as follows: In section 2, we define the kinematics of the process B — * 
X s £ + £~ and introduce the quantities of dynamical interest in the framework of an effective Hamil- 
tonian. Section 3 describes the wave-function effects in the FM model in the hadron energy and 
hadronic invariant mass spectra. The resulting inclusive hadron spectra in B — » X s £ + l~ , including 
the long-distance effects in terms of the J/ip, ip', ... resonances, are presented. Comparison with the 
J/^-momentum spectrum measured by the CLEO collaboration is shown, constraining the FM model 
parameters. Theoretical uncertainties in the effective coefficients Re Cg S and |C| ff | from the various 
prescriptions are displayed here and the resulting spectra in hadron energy, dilepton invariant mass 
and the FB-asymmetry are presented. Hadronic spectral moments are calculated in the FM model 
taking into account the LD contributions. Comparison with the corresponding quantities derived from 
the SD-contribution alone using the HQET and FM models in fiJ0| is also presented. In section 4, the 
effects of the experimental cuts used in the CLEO analysis of B — ► X s £ + £~are studied and the resulting 
spectra are presented in terms of several figures. Estimates of the branching ratios B(B — ► X s £ + £~) 
for £ = e, fx are also presented here, together with estimates of the survival probability for the CLEO 
cuts, using the FM model. Section 5 contains a summary of our work and some concluding remarks. 

2 The Decay B — > X s £ + £ in the Effective Hamiltonian Approach 
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2.1 Kinematics 

The kinematics for the decay in question at the partonic and hadronic level are defined in Hence, 
we will be short here. The parton level kinematics is given by 

b(p b ) -> s( Ps )(+g( Pg )) + £ + ( P+ ) + r (p_) , (1) 

where g denotes a gluon from the 0(a s ) correction. The corresponding kinematics at the hadron level 
is defined as: 

B(p B ) -> X S ( PH ) + £+(p + ) + r (p_) . (2) 

We define the momentum transfer to the lepton pair and the invariant mass of the dilepton system, 
respectively, as 

q = p++p- , (3) 
s = q 2 . (4) 

Further, we define a 4-vector v, which denotes the velocity of both the 6-quark and the 5-meson, 
Pb = mtfv and pb = tubv. The hadronic invariant mass is denoted by Sh = v\ an d Ejj denotes the 
hadron energy in the final state. The corresponding quantities at parton level are the invariant mass so 
and the scaled parton energy xq = ^l. In parton model without gluon bremsstrahlung, this simplifies 
to so = vn s and xq becomes directly related to the dilepton invariant mass xq = l/2(l-s + 7r^). Here 
and in what follows, the dimensionless variables with a hat are related to the dimensionful variables 
by the scale m&, the 6-quark mass, e.g., s = rh s = ^ etc. From momentum conservation the 
following equalities hold in the 6-quark, equivalently 5-meson, rest frame (v = (1,0,0,0)): 

x = 1-vq, 

s = 1 - 2v ■ q + s , (5) 

Eh = rriB — v ■ q , 

Sh = m B — 2m B v ■ q + s . (6) 

The relation between the kinematic variables of the parton model and the hadronic states can be seen 
in !!■ 

2.2 Matrix element for the decay B — > X s t r t~ 

The effective Hamiltonian obtained by integrating out the top quark and the bosons is given as 

6 

Ci{pL)Oi + Cririj^Eaa^mbR + m s L)b a F^ 

i=l 

^ ^ 1 

(7) 



H ef f(b - s + X, X = 7, £+£-) = -^-V t * s V tb 



where L and R denote chiral projections, L(R) = 1/2(1 =p 75), V« are the CKM matrix elements and 
the CKM unitarity has been used in factoring out the product V t *Vtb- The operator basis is taken 
from [|l(]] , where also the Four-Fermi operators 0\ , . . . , Og and the chromomagnetic operator Og can 
be seen. Note that Og does not contribute to the decay B — > X s £ + £~in the approximation which 
we use here. The Cj(//) are the Wilson coefficients, which depend, in general, on the renormalization 
scale n, except for Cio- 

The matrix element for the decay B — > X s £ + £~can be factorized into a leptonic and a hadronic 
part as 

M(B^X s £+£-) = ^-V? s V tb (T\L L » + T R ,L R ») , (8) 

with 



T L/R 



(cf(s) T Cio + 2Cf j) + 2m s cf^ j-L 



b. 



(9) 
(10) 



rr 

The effective Wilson coefficient Cg (s) receives contributions from various pieces. Since the resonant 
cc states also contribute to Cg^(s), the contribution given below is just the perturbative part: 



cf(s)\ pert = C 9V (s)+Y(s). 



(11) 



Here rj(s) and Y(s) represent the 0(a s ) correction [21| and the one loop matrix element of the Four- 
Fermi operators p2]j2^ ], respectively. While Cg is a renormalization scheme-dependent quantity, this 
dependence cancels out with the corresponding one in the function Y(s) (the value of £, see below). 
To be self-contained, we list the two functions in Cg (s): 



Y(S) = g(m c ,s)(3C l + C 2 + SC 3 + C 4 + 3C 5 + C 6 ) 

~g(l, s) (4 C 3 + 4 C A + 3 C 5 + C 6 ) - ±g(0, s) (C 3 + 3 C 4 ) 

+1 (3C 3 + C 4 + 3C 5 + C 6 ) - ^ (3Ci + C 2 - C 3 - 3C 4 ) , 
17(5) = l + ^W), 



7T 



(NDR), 
-1 (HV), 



(12) 
(13) 

(14) 



8 m b 

9 ~~ 

G(l-y)(ln 



8 4 
mz H 1 — y 

27 r 
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(2 + yW|l-y 



1 + yT^ 

i-yr^ 



nr) + Q(y — l)2arctan 



rn « 8 8 m fe 

*(<>,*) = 57 -9 My) 



4, A 4. 
— ins H — itt , 
9 9 ' 



(15) 
(16) 



where y = 4z /§, and 
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4 

3 L «2(s) 



— In s ln(l 



5 ± 4s 
3(1 + 2s 



- ln(l 



2s(l±s)(l-2s) 1 „ 5 + 9s 
Ins + 



6s 2 



3(1 



£^2 



(1 + 2s) 



6(l-s)(l±2s) 



(17) 



Above, (NDR) and (HV) correspond to the naive dimensional regularization and the 't Hooft-Veltman 
schemes, respectively. Note that the function g(m c , s)(3C\ + C<i ± ...) given above is the perturbative 
contribution to the effective coefficient Cg (s) from the cc loop, to which we have referred in the 
introduction and to whose discussion we shall return in section 3. The other Wilson coefficients in 



leading logarithmic approximation can be seen in [22 



Parameter 


Value 


mw 


80.26 (GeV) 


m z 


91.19 (GeV) 


sin 2 9w 


0.2325 


m s 


0.2 (GeV) 


m c 


1.4 (GeV) 


m b 


4.8 (GeV) 


m t 


175 ± 5 (GeV) 




+m h 


m »-m b b /2 


A (5) 
ly QCD 


0.214i°;«* (GeV) 


a- 1 


129 


a, s (m z ) 


0.117 ±0.005 


B s i 


(10.4 ±0.4) % 


Ai 


-0.20 (GeV 2 ) 


A 2 


±0.12 (GeV 2 ) 



Table 1: Default values of the input parameters and errors used in the numerical calculations. 



Ci 


c 2 


c 3 


c 4 


c 5 


c 6 




c 9 


Cio 


c<f>) 


-0.240 


±1.103 


±0.011 


-0.025 


±0.007 


-0.030 


-0.311 


±4.153 


-4.546 


±0.381 



Table 2: Values of the Wilson coefficients used in the numerical calculations corresponding to the 
central values of the parameters given in Table |^. For Cg we use the NDR scheme. 



3 Hadron Spectra in B — > X s £ + £ in the Fermi Motion Model Includ- 
ing Long-Distance Contribution 

In this section, we study non-perturbative effects associated with the bound state nature of the B 
hadron and of the cc resonances B — > X s ± (J/ift, t//, ■■■) — ► X s l + i~ (the LD contribution) on the 



hadronic invariant mass and hadron energy distributions in the decay B — ► X s £ + £ . Wave-function 
effects are studied in the FM model Q and for the resonant part we use data. The FM model 



parameters are then constrained from the measured J/^-momentum spectrum [17]. Sensitivity of 
the inclusive spectra in B — > X s £ + £~on the assumed resonant and the perturbative cc-contribution 
is presented. The hadronic spectral moments (Xjj), with X = S,E and n = 1,2, are calculated 
numerically in the FM model including the LD-effects. For the sake of comparison, the corresponding 
quantities calculated for the SD-contribution using HQET and the FM model from M are also given. 

3.1 Hadron spectra in B — > X s £ + £~\\\ the Fermi motion model 

The FM model often invoked in phenomenological studies of B decays [|| has two parameters pf, the 
Fermi momentum of the 6-quark, and the spectator quark mass m q . Energy- momentum conservation 
requires the 6-quark mass to be a momentum-dependent parameter determined by the constraint: 



m b (p)=m B + m q - 2m B ^p 2 + m 2 ; p = \p\ . (18) 

The 6-quark momentum p is assumed to have a Gaussian distribution, denoted by 4>(p), which is 
determined by pp 

4 — p 2 

Hp) = —f= — ? exp( — -) , (19) 

^-Kpp A PF Z 

with the normalization / °° dpp 2 4>(p) = 1. In this model, the HQET parameters A and Ai, repre- 
senting, respectively, the binding energy and the kinetic energy of the 6-quark inside a B meson, are 
calculable in terms of pp and m q with 



A = J dpp 2 4>{p)\jm\ + p 2 , 



3 
2 1 



Ai = - / dppU(p) = -^V 2 f- (20) 



In addition, for m q = 0, one can show that A = 2pi?/y / 7r. There is, however, no analog of A2 in the 
FM model. For subsequent use in working out the normalization (decay widths) in the FM model, we 
also define an effective 6-quark mass by 



m b — 







dpp 2 m b {pf(t){p)) 1 ^ . (21) 

With the quantities , Ai and A defined above, the relation 

m B = mf + A - A!/(2m^ ff ) , (22) 

is found to be satisfied in the FM model to a high accuracy (within 0.7%), which is shown in Table 
H for some representative values of the HQET parameters and their FM model equivalents. We shall 
use the HQET parameters A and Ai to characterize also the FM model parameters, with the relations 
given m eqs. @ and @ and in Table |. 
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/"IV /T "V T IV T ... "\ 7"\ 

PF,m q (MeVjMeV j 


eft 1 ' r\ t r\ 

(GeV) 


Ai (GeV ) 


A (GeV) 


(450,0) 


4.76 


-0.304 


0.507 


(252,300) 


4.85 


-0.095 


0.422 


(310,0) 


4.92 


-0.144 


0.350 


(450, 150) 


4.73 


-0.304 


0.534 


(500, 150) 


4.68 


-0.375 


0.588 


(570, 150) 


4.60 


-0.487 


0.664 



eff 

Table 3: Values of non perturbative parameters m b , Ai and A defined in the FM model for different 
sets of the FM model parameters (pF,m q ). 



3.2 Long-distance resonant contribution in B — > X s £ + £ 

If not stated otherwise, we shall follow the procedure adopted in fTofl , in which the long-distance 
(LD) resonance effects in the decay B — * X s £ + £~ (specified below) are added with the perturbative cc 
contribution expressed through the function g(m c ,s) in section 2 (see, eq. (pi])). Thus, 

Cf(s) = C 9V (s) + Y(s) + Y res (s) . (23) 

The function Y res (s) accounts for the charmonium resonance contribution via B — > X s (J/^f, \&', ...)—► 

TTT(Vi 



X s £ + £ for which we take the representation [11|, 

3 



Y r es (s) 



E 



e + £-)M v . 

M v *-sm b 2 -iM v T Vi 



(24) 



Vi=iP(ls),...,^(6s) 

where = 3Ci +C 2 + 3C 3 + C 4 + 3C 5 + C 6 . We adopt k = 2.3 for the numerical calculations @. 
This is a fair representation of present data in the factorization approach |25]] ; also the phase of k, 
which is fixed in eq. fl24|), is now supported by data which finds it close to its perturbative value ||. 
Note that in this approach, the effective coefficient Cg^(s) has a s-dependence, which is not entirely 
due to the propagators in the function Y res (s) as also the perturbative cc contribution g{rn c ,s) is a 
function of s. In the resonant region, the perturbative part is not noticeable due to the fact that the 
resonant part in Cg^(s) completely dominates. However, when the cc pair is sufficiently off-shell, the 
s-dependence of the function Cg^(s) is not (and should not be) entirely determined by the cc resonant 
contribution. This is the motivation of the representation in eq. (23). We shall later evaluate the 
uncertainties in various distributions arising from varying definitions of the perturbative contribution 
g(m c , s) as well as the precise form of the resonating contribution. 

3.3 Constraints on the FM model parameters from existing data 

The FM model parameters pf and m q (equivalently Ai and A) for the SD-contribution can, in principle, 
be determined from an analysis of the energy spectra in the decays B — * X u £f£ and B — > X s + 7, as all 

7 




Figure 1: Momentum distribution of J/i/j in the decay B — ► X S J /ip in the FM model. The solid, 
dotted, dashed curve corresponds to the parameters (Ax, A) = (—0.3, 0.5), (—0.3, 0.53), (—0.38, 0.59) in 
(GeV 2 , GeV), respectively. The data points are from the CLEO measurements W 



of them involve the decay of a b quark into (an almost) massless (u or s) quark. However, the quality 
of the present data on B — > X s + 7 does not allow to draw very quantitative conclusions, and hence 
we vary the parameters in a reasonable range giving a satisfactory description of the photon energy 
spectrum and show the resulting uncertainties on the hadron spectra in B — > X s i + l~ . For estimating 
the spectra from the LD contribution involving the transition b — ► c, the parameters of the FM model 
can be constrained from the lepton energy spectrum in the decay B — > X c tv^ and from the shape 
of the J/tp- and tp'- momentum distributions in the decays B — * X s (J/ip,ip'). We review below the 
presently available analyses of the photon- and lepton-energy spectra in B decays in the FM model 
and also present an analysis of the J/^-momentum spectrum in B — > X s J/ip. 

• Analysis of the photon energy spectrum in B — » X s + 7 

The photon energy- and invariant hadronic mass distributions in B — > X s ^f were calculated in the 



FM model using the leading order (in a s ) corrections in ref. [26]. These spectra were used in the 
analysis of the CLEO data on B — > X s + 7 |27|, in which the values pf = 270 ± 40 MeV suggested 
by the analysis of the CLEO data on B — > X£u£ were used, together with the effective 6-quark mass 

rr 

= 4.87 db 0.10 GeV, which gave reasonable fits of the data. We translate these parameters in 
terms of Ai and A using the relations given in eqs. ( ]20j ) and (f2^), yielding 

Ai = -0.11+g;g^ GeV 2 , A = 0.40 ± 0.1 GeV . (25) 



10 7 g- 
10 6 r 

k 10 5 =- 

= 

O 1Q4 r 

ljj : 
\ 1000 - 

1 100 r 



A 10 r 

m c 



1 r 

0.1 - 





i i L 

2 

E H [GeV] 



Figure 2: Hadron energy spectrum in B — > X s £ + £~ including the resonant (J/ip,ip' , ...) anrf perturba- 
tive contributions in the Fermi motion model (dotted curve) for (Ai, A) = (—0.1 GeV 2 , OA GeV), and 
in the parton model (long-short dashed curve) for mi, = 4.85 GeV. 

The same data was fitted in ref. in the FM model, yielding (pp,m q ) = (0.45 GeV,0 GeV) as the 
best-fit solution, with (pF,m q ) = (0.310 GeV, 0.3 GeV) differing from the best-fit solution by one unit 
in x 2 - The best-fit values translate into 



\ l = -0.3 GeV 2 , A = 0.5 GeV . 



(26) 



Within the indicated errors, the values given in eqs. (|25| ) and (^g) are compatible. 

• Analysis of the lepton energy spectrum in B — > Xiv^ 

A fit of the lepton energy spectrum in the semileptonic decay B — » Xlvi in the context of HQET has 
been performed in ref. Using the CLEO data |2q| , the authors of ref. || find: 



A I 



-0.19 ± 0.10 GeV , A = 0.39 ± 0.11 GeV . 



(27) 



Since the FM model and HQET yield very similar lepton energy spectra (apart from the end-point), 
one can take the analysis of || also holding approximately for the FM model. 

• Analysis of the J/^-momentum spectrum in B — > X s J/ip 

An analysis of the J/'i/'-momentum spectrum in B — > X s (J/ijj, ip') measured by the CLEO collaboration 
[jnj] has been reported in ref. [18] using the FM model. The authors of ref. ]l8[ addressed both the 
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Figure 3: Hadronic invariant mass spectrum in B — > X s £ + £ including the perturbative and resonant 
(J/ifj,if)' , ...) contributions in the Fermi motion model. The solid, dotted and dashed curves correspond 
to the parameters (Ai,A) = (-0.3, 0.5), (-0.1, 0.4), (-0.15, 0.35) in (GeV 2 , GeV), respectively. 



shape and normalization of the J/^-data, using the non-relativistic QCD (NRQCD) formalism for 
the inclusive color singlet and color octet charmonium production in B — > X s J/ijj and the FM model 
[||. The preferred FM parameters from this analysis are: (pp,m q ) = (0.57 GeV, 0.15 GeV), where 
rriq only plays a role in determining the position of the peak but otherwise does not influence the 
small momentum tail of the J ftp momentum distribution. This yields values of the parameter pp 
which are consistent with the ones obtained in ref. [29] pp = 0.54^Q'Jg, GeV based on an analysis of 
the CLEO data on B — > X£v£ p8[| . The central values of pp in |29] as well as in |l8| correspond to 
m en ^ ^ q GeV, which is on the lower side of the present theoretical estimate of the m& pole mass, 
namely m b = 4.8 ± 0.2 GeV 

We have redone an analysis of the J/ -^-momentum distribution in the FM model which is shown 
in Fig. [I]. As also discussed in [18], the low-momentum J/ip, in particular in the region < 0.6 

GeV, are problematic for inclusive decay models, including also the FM model (see Fig. 1). The 
measured -spectrum appears to have a secondary bump; an inclusive spectrum behaving as a 

Gaussian tail or having a power-like behavior oc \kj/^\~ 5 in this region is hard put to explain this 
data. There are also suggestions in literature [^TJ] that the spectrum in this region is dominated by the 
three-body decay B — * J /ipKp and hence the bump reflects the underlying dynamics of this exclusive 
decay. In view of this, we have taken out the first six points in the \ow-\kj/^\ spectrum and fitted the 
FM model parameters in the rest of the \kj/^ [-spectrum. The three curves shown correspond to the 
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rr 

Figure 4: The real part (a) and the absolute value (b) of C$ (s) are shown as a function of s, where 

rr 

C| n (s) = C 9 rt(s)+Y(s)+Y res (s). The solid curve corresponds to Y(s) calculated using the perturbative 



cc contribution g(rh c , s) given in eq. (15), and the dotted curve corresponds to using g(rh c , s) in eq. 



with R res (s) calculated in both cases using eq. (24). The dashed curve corresponds to using the Kriiger- 
Sehgal approach fTl] discussed in the text. 



FM model parameters (pF,m q ) = (0.45 GeV,0 GeV) (solid curve), (pF,m g ) = (0.45 GeV, 0.15 GeV) 
(dotted curve) and (j>F,m q ) = (0.50 GeV, 0.15 GeV) (dashed curve). They all have reasonable x 2 ■, 
with x 2 /dof = 1.6, 1.6 and 1.1, respectively. Excluding also the seventh lowest point, the x 2 improves 
marginally, with the resulting x 2 being x 2 /dof = 1.4,1.4 and 0.94. Including the sixth point, the 
fits become slightly worse. However, they are all acceptable fits. It is interesting that the best-fit 
solution of the photon energy spectrum in B — ► X s + 7, (pp, m q ) = (0.45 GeV, GeV) | |T9[ | , is also an 
acceptable fit of the |/cj/^,|-data. The corresponding Ax, A and values from these two analyses are 
compatible within drier with the HQET-based constraints from the semileptonic B decays ]|], quoted 



above. Thus, the values in eq. (26) appear to be a reasonable guess of the FM model parameters. 
But, more importantly for the present study, the phenomenological profile of the LD contribution 
B — > X s (J/tp, ip' , ...) — > X s £ + £~ presented here is certainly consistent with present data and theoretical 
constraints. 

3.4 Effects of the Lorentz boost on the hadron spectra in B — > X s £ + £~ 

We now discuss the i?-meson wave function effects in the FM model on the hadron spectra in B — > 
X s £ + £~ . Since the resonances in B — > X s £ + £~&ie in the dilepton invariant mass variable s and not 
in Sh, and noting that neither Eq (partonic energy) nor Eh are Lorentz-invariant parameters, it 
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Figure 5: Hadron energy spectrum in B —* X s £ + £~ including the resonant and perturbative con- 
tributions in the Fermi motion model. In (a), the FM model parameters are fixed at (Ai,A) = 
(—0.1 GeV 2 , OA GeV). The almost overlapping curves differ in the parametrization of Cg (s) as 
depicted in Fig. (4). The solid curve is obtained using eq. for g(fh c ,s), the dotted curve is based 
on g(rh c , s) given in eq. ftBoi), with R res (s) calculated in both cases using eq. $24V, and the dashed curve 
corresponds to the Kriiger-Sehgal approach ^|/. In (b), the solid, dotted and dashed curves correspond 
to the parameters (Ai,A) = (-0.3, 0.5), (-0.1, 0.4), (-0.15, 0.35) in (GeV 2 , GeV), respectively. 

is expected on general grounds that the effect of the Lorentz boost in the FM model on Eh- and 
<Sf/-distributions will be more marked than what was found on the invariant dilepton mass spectrum 
in [ 10 1 . We recall that for the dilepton invariant mass, the Lorentz boost involved in the FM model 
(Doppler shift) leaves the spectrum invariant and there is only a residual effect due to the fact that the 
6-quark mass in the parton model nif, and the effective 6-quark mass in the FM model (called W(p) 
in |l(| and here), are different quantities. This difference (W(p) — m&) (mass defect) smears the 
dilepton invariant mass distribution, but being a subleading effect in l/m& this effect is small. Not so 
in the hadron spectra. In the hadron energy spectrum, the cc-resonances, which are narrowly peaked 
in the parton model, are broadened by the Lorentz boost of the FM model. To show this, the hadron 
energy spectrum in the FM model (dotted curve) is compared with the spectrum in the parton model 
(long-short dashed curve) in Fig. ^ for identical values of mj, and mf^, taken as 4.85 GeV. In terms 
of the hadronic invariant mass, one finds that the resonant structure is greatly smeared. The reason 
for this behavior is that each q 2 -bin contributes to a range of Eh and Sh- The different q ,2 -regions 
overlap in Sh resulting in a smearing of the resonances over a wide range. This can be seen in Fig. |3] 
for the hadronic invariant mass. Various curves illustrate the sensitivity of this spectrum on the FM 
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(a) 



(b) 




5 10 15 

q 2 [GeV 2 ] 




10 15 

q 2 [GeV 2 ] 



Figure 6: Dilepton invariant spectrum (a) and the (normalized) Forward- Backward asymmetry (b) in 
B — > X s l + £~ including the resonant and perturbative contributions in the Fermi motion model. The 
FM model parameters are fixed at (Ai, A) = (—0.1 GeV 2 , OA GeV). The solid curve is obtained using 
eq. (j73p for g(rh c ,s), the dotted curve is based on g(m c ,s) given in eq. ft23j), with R res (s) calculated in 



both cases using eq. (]^4j), and the dashed curve corresponds to the Kriiger-Sehgal approach 



model parameters. 

3.5 Ambiguities in adding LD and SD contributions in B — > XJ + i~ 

Since we are simply adding the short-distance (SD) and resonant charmonium amplitudes, it can not 

rr 

be ruled out that possibly some double counting has crept in in the coefficient Cf (s) (once as a 
continuum cc contribution and then again as J/ip,ip', ... resonances). The question is whether the 



addition of the cc-continuum and resonating pieces as being done here and in [1C] compromises the 
resulting theoretical precision significantly. This can only be studied by comparing the theoretical 
scenario in question with other trial constructions advocated in the literature. For example, one 
could retain in the perturbative function g(rh c , s) just the constant part in s by replacing g(rh c , s) by 
g(rh c , s), where 



9 V n ' 



— In m r H 

9 27 



(28) 



This function (with fi = m&) has been proposed in [|l4|,|l^] as an alternative representation of the 
cc perturbative contribution and represents the (minimal) short-distance contribution. To study the 
difference numerically, we plot both the real part ReCg^(s) and the absolute value |Cg (s)| as func- 
tions of s in Fig. |3| by using the complete perturbative expression for g(rh c , s) in eq. fli~2| ) and g(rh c , s) 
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given in eq. (pq). In both cases, the resonant contributions are included using eq. (^4|). As a third 
parametrization of Cg (I), we use the the approach of Kriiger and Sehgal | |14| |, based on dispersion 
relation Q : 

Im g(m c ,s) = ^R^s) , 

Re g(m e , s) = -§ lnm c - - + | / % da' . (29) 

9 9 3 7 4 m2 s'(s' - s) 

The cross-section ratio R^^{s) in this approach is expressed as 

R^) = RcUs) + Rfes(s), (30) 

where Rcont(s) an d -^res(s) denote the contribution from the continuum and the narrow resonances, 
respectively. For the narrow resonances, the Breit-Wigner form given below in eq. ( |3ll) is used, whereas 
for the continuum part a parametrization of the e + e~ annihilation data, taken from |H]], is used. 

^res^J - K 2^ 7^^~2 ™2 N2 , m 2 t^2 ' ^ 



total 



Note that the authors of [14| use k = 2.35 in their numerical analyses, which is slightly different than 
the one used in [10]. More importantly, this parametrization of R res (s) has a different s-dependence 
than the one given by eq. (^) and it incorporates the off-shell dependence of the effective 7* — ► V 



vertex, discussed in fll3|. 

A number of comments on the curves shown in Fig. || and the resulting hadron energy and hadronic 
invariant mass distributions is in order. 

• The results for both Re Cg^(s) and the absolute value |Cg^(s)| plotted as functions of s show 
that the functions corresponding to [|l4|] are lower than the ones used in which, in turn, are 
lower than the ones obtained with the prescription given in [15]. However, one sees from Figs. |] 
(a) and || (b) that the differences between these functions are numerically small. 

• The three parametrizations discussed above give almost identical hadron spectra. The differences 
between these approaches in the .E#-spectrum are already difficult to see, as shown in Fig. || (a) 
where we have plotted the iij^-spectra in the FM model for the three parametrizations; the effect 
on the hadronic invariant mass is even less noticeable and hence is not shown. Quantitatively, the 
maximum difference in the hadron energy and the hadronic invariant mass spectra is 12.1(4.5)% 
and 4.1(2.5)%, respectively. The difference between the approaches in |l(| and [14| is larger, 
as given by the first numbers, than the one between Jlc| ] and [15], given by the numbers in 
parentheses. However, these maximum differences occur only over a rather limited part of the 
phase space. 

1 The lower limit of the integral in eq. ( ps| ) must be below rh 2 j^ to include the dominant lowest cc resonance contri- 
butions. We use 4mJ in consultation with F. Kriiger. 
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• Other uncertainties on the hadronic distribution are much larger, see, for example, Fig. |5| (b) 
showing the sensitivity of the hadron energy spectra on the 5-meson wave-function parameters. 
In comparison, the cc-resonance/continuum related ambiguity is numerically small. 

The results presented in Figs. ^ and || are the principal phenomenological results derived by 
us in this paper for the inclusive hadronic invariant mass and hadron energy spectra in the decay 
B — > X s £ + £~ , respectively, and are of direct experimental interest. The £f/-distribution and moments 
depend on the FM model parameters pf and m q . In the HQET approach, they depend on the 
parameters A and Ax- Since these parameters are already constrained by present data, the decay 
B —* X s £ + £~can be gainfully used to determine them more precisely. This was discussed in [[!],§]. 
Here, we have estimated the influence of the LD-resonant contribution. 

3.6 Numerical Estimates of the Hadronic Moments in FM model and HQET 

The similarity of the first two hadronic moments (Xjj), with X = E, S and n = 1,2, involving 
the SD-contribution in the decay B —* X s £ + £~'va. the HQET and FM model descriptions was shown 
quantitatively in 0. We include these numerical results here for comparison with the corresponding 
moments calculated including the LD-contribution in the FM model using the spectra which we have 
presented above. The moments based on the SD-contribution are defined as: 

(X n H ) = ( / X n H ^-dX H )jB for X = S, E . (32) 

The moments (X|^) gc are defined by taking into account in addition to the SD-contribution also the 
contributions from the cc resonances. The values of the moments in both the HQET approach and 
the FM for n = 1,2 are shown in Table |], with the numbers in the parentheses corresponding to 
the former. They are based on using the central values of the parameters given in Table |l], except 
for the values of Ai and A which are explicitly stated. The correspondence between the FM model 
and HQET parameters is given in eqs. (|20|). As already stated in Q, both the HQET and the FM 
model lead to strikingly similar results for the SD-contribution based hadronic moments shown in this 
table. However, the moments {X^) Sc with X = S,E are significantly lower than their SD-counterparts 
{Xft) calculated for the same values of the FM model parameters. This shows that the cc resonances 
are important also in moments. The hadronic invariant mass spectra in B — > X s £ + £~foi both the 
SD and inclusive (LD+SD) contributions are expected to be dominated by multi-body states, with 
(S H ) - (1-5 - 2.1) GeV 2 and (S H )cc (1.2 - 1.5) GeV 2 . 
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(Sh) 


(Sh)cc 


(S%) 


{S"h)cc 


(Ai, A) in (GeV 2 , GeV) 


(GeV 


2 ) 


(GeV 


4 ) 


(-0.3,0.5) 


2.03 (2.09) 


1.51 


6.43 (6.93) 


3.10 


(-0.1,0.4) 


1.75 (1.80) 


1.36 


4.04 (4.38) 


2.17 


(—0.14, 0.35) 


1.54 (1.49) 


i 1 n 

1.19 


3.65 (3.64) 


1.92 




(Eh) 


(Eh)cc 


(E 2 H ) 


/ TT'2 \ 

\E H )cc 


(Ai, A) in (GeV 2 , GeV) 


(GeV) 


(GeV 


2 ) 


(-0.3,0.5) 


2.23 (2.28) 


1.87 


5.27 (5.46) 


3.52 


(-0.1,0.4) 


2.21 (2.22) 


1.85 


5.19 (5.23) 


3.43 


(-0.14,0.35) 


2.15 (2.18) 


1.84 


4.94 (5.04) 


3.39 



Table 4: Hadronic spectral moments for B — > X s /j, + /j, in the Fermi motion model (HQET) for the 
indicated values of the parameters (Ai,A). 

4 Branching Ratios and Hadron Spectra in B — > X s £ + £ with Cuts 
on Invariant Masses 

In experimental searches for the decay B —* X s £ + £~, the short-distance contribution (electroweak 
penguins and boxes) is expected to be visible away from the resonances. So, cuts on the invariant 
dilepton mass are imposed to stay away from the dilpeton mass range where the charmonium reso- 
nances J/ip and ij)' are dominant. For example, the cuts imposed in the recent CLEO analysis |@] 
given below are typical: 

cut A : q 2 < (m J/4) - 0.1 GeV) 2 = 8.98 GeV 2 , 
cut B : q 2 < (mj/^ - 0.3 GeV) 2 = 7.82 GeV 2 , 

cut C : q 2 > (mp + 0.1 GeV) 2 = 14.33 GeV 2 . (33) 

The cuts A and B have been chosen to take into account the QED radiative corrections as these effects 
are different in the e + e _ and modes. In the following, we compare the hadron spectra with and 

without the resonances after imposing these experimental cuts. For the low-g 2 cut for muons (cut ^4), 
the hadron energy spectra and the hadronic invariant mass spectra are shown in Fig. [?] (a), (b) and 
Fig. |8| (a), (b), respectively. The results for the low-g 2 cut for electrons (cut B), are shown in Fig. |7| 
(c), (d) and Fig. || (c), (d), respectively. Finally, the hadronic spectra for the high-q 2 cut (cut C) for 
e + e~ and fi + can be seen in Fig. [?] (e), (f) for the hadronic energy and in Fig. || (e), (f) for the 
hadronic invariant mass. We see that the above cuts in q 2 greatly reduce the resonance contributions. 
Hence, the resulting distributions essentially test the non-resonant cc and short-distance contributions. 
These figures will be used later to quantify the model dependence of the integrated branching ratios 
in B X s £+e~. 
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As mentioned in [J/J, the dominant BB background to the decay B — > X s ^ + £~comes from two 
semileptonic decays of B or D mesons, which produce the lepton pair with two undetected neutrinos. 
To suppress this BB background, it is required that the invariant mass of the final hadronic state 
is less than t = 1.8 GeV, which approximately equals mo. We define the survival probability of the 
B — > X s ^ + ^ _ signal after the hadronic invariant mass cut: 

S(t) = ( [\ ^-dS H )/B , (34) 

Jm 2 x CLOH 

and present S(t = 1.8 GeV) as the fraction of the branching ratio for B — > X s £ + £~ surviving these cuts 
in Table ||. To estimate the model dependence of this probability, we vary the FM model parameters. 
Concentrating on the SD piece, we note that the effect of this cut alone is that between 83% to 92% 
of the signal for B — ► X s fi + ^~ and between 79% to 90% of the signal in B — > X s e + e~ survives, 
depending on the FM model parameters. The corresponding numbers for the inclusive spectrum 
including the SD and LD contribution is 96% to 99.7% for both the dimuon and dielectron case. 
This shows that while this cut removes a good fraction of the BB background, it allows a very large 
fraction of the B — > X s .£ + i!~ signal to survive. However, this cut does not discriminate between the 
SD and (SD+LD) contributions, for which the cuts A - C are effective. The numbers for the survival 
probability S(t = 1.8 GeV) reflect that the hadronic invariant mass distribution of the LD-contribution 
is more steep than the one from the SD contribution. 

With the additional cut A (B) imposed on the dimuon (dielectron) invariant mass, between 57% 
to 65% (57% to 68%) of the B — ► X s £ + £~ signal survives the additional cut on the hadronic invariant 
mass for the SD contribution. However, as expected, the cuts A and B result in drastic reduction 
of the inclusive branching ratio for the decay B — ► X s £ + £~, as they effectively remove the dominant 
cc-resonant part. In this case only 0.8% to 0.9% (1.0% to 1.2%) of the inclusive signal survives for 
the cut A (B). The theoretical branching ratios for both the dielectron and dimuon cases, calculated 
using the central values in Table [l] and the indicated values of Ai and A are also given in Table ||. As 
estimated in [|J, the uncertainty on the branching ratios resulting from the errors on the parameters 
in Table [l] is about ±23% (for the dielectron mode) and ±16% (for the dimuon case). The wave- 
function-related uncertainty in the branching ratios is smaller, as can be seen in Table [|. This gives 
a fair estimate of the theoretical uncertainties on the partially integrated branching ratios from the 
S-meson wave function and cc resonances. With the help of the theoretical branching ratio and the 
survival probability S(t = 1.8 GeV), calculated for three sets of the FM parameters, the cross section 
can be calculated for all six cases: 

(i) no cut on the dimuon invariant mass [(SD) and (SD + LD)], (ii) no cut on the dielectron invariant 
mass [(SD) and (SD + LD)], (iii) cut A on the dimuon invariant mass [(SD) and (SD + LD)], (iv) cut B 
on the dielectron invariant mass [(SD) and (SD + LD)], (v) cut C on the dimuon invariant mass [(SD) 
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and (SD + LD)], (vi) cut C on the dielectron invariant mass [(SD) and (SD + LD)]. This table shows 
that with 10 BB events, O(70) dimuon and 0(100) dielectron signal events from B — » X s ^ + £~should 
survive the CLEO cuts A and B, respectively, with m(X s ) < 1.8 GeV. With the cut C, one expects 
an order of magnitude less events, making this region interesting for the LHC experiments which will 
have much higher BB statistics. Given enough data, one can compare the experimental distributions 
in B — ► X s £ + £~ directly with the ones presented here. The phenomenological success of the FM model 
in describing the energy spectra in B decays and its close proximity to HQET make us confident that 
the hadron spectra in B — » X s ^ + £~presented here should be good descriptions of the data. 

4.1 Hadronic Spectral Moments with Cuts in the FM 

We have calculated the first two moments of the hadronic invariant mass in the FM model by imposing 
a cut Sjj < t 2 with t = 1.8 GeV and an optional cut on q 2 . 

<ss) - </I s »w S? dsrf > for - = ^ 2 ■ < 35 » 

Here the subscript cutX indicates whether we evaluated (Sh) and (S H ) with the cuts on the invariant 
dilepton mass as defined in eq. (|33|), or without any cut on the dilepton mass. The results are collected 
m Table § The moments given in Table ^ can be compared directly with the data to extract the 
FM model parameters. The entries in this table give a fairly good idea of what the effects of the 
experimental cuts on the corresponding moments in HQET will be, as the FM and HQET yield very 
similar moments for equivalent values of the parameters. The functional dependence of the hadronic 
moments on the HQET parameters taking into account the experimental cuts still remains to be 
worked out. 

In the last row of Table || the value in percentage refers to the maximum uncertainty in (Sjj), with 
n = 1, 2, resulting from different approaches to include the resonant cc effects. We have calculated (Sh) 
and (Sjj) with all the cuts mentioned above for the three approaches for fixed (Ai,A) = (—0.1,0.4) 
in GeV 2 , GeV. Thus, for example, (S H ) = (1-77 ± 0.90%) GeV 2 for Cut A. This uncertainty is much 
below the one due to the variations in the parameters Ai and A. Hence, the measurement of the 
hadronic moments can be used to determine these parameters. 

5 Summary and Concluding Remarks 

We summarize our results: 

• We have presented the hadron spectra and moments in B — » X s £ + £~ including the resonant- 
cc contribution in the Fermi motion model Q. This complements the description of the final 
states in B — > X s £ + £~ presented in |i~0[| , where the dilepton invariant mass spectrum and FB 
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FM parameters 


B ■ 10" 6 


B ■ 10- 6 


No s-cut 


No s-cut 


cut A 


cut B 


cut C 


cut C 


(Ai, A) in (GeV 2 , GeV) 


1 


e^e 


1 




1 


e^e 


1 


e^e 


(-0.3,0.5) [SD] 


5.8 


8.6 


83% 


79 % 


57% 


57% 


6.4% 


4.5% 


(-0.1,0.4) [SD] 


5.7 


8.4 


93% 


91 % 


63% 


68% 


8.3% 


5.8% 


(-0.14,0.35) [SD] 


5.6 


8.3 


92% 


90 % 


65% 


67% 


7.9% 


5.5% 


(-0.3,0.5) [SD + LD] 


562.5 


563.9 


96% 


96 % 


0.8% 


1.0% 


0.06% 


0.06% 


(-0.1,0.4) [SD + LD] 


564.0 


565.6 


99.7% 


99.7% 


0.8% 


1.1% 


0.08% 


0.08% 


(-0.14,0.35) [SD + LD] 


566.5 


568.2 


99% 


99 % 


0.9% 


1.2% 


0.08% 


0.08% 



Table 5: Branching ratios and survival probabilities for B — > X s l + £~ , £ = //,e for different FM 
model parameters evaluated from the SD and [SD + LD] contributions. The branching ratios without 
experimental cuts are given in the second and third columns. The values given in percentage in the 



fourth to eleventh columns represent the survival probability S(t = 1.8 GeV) defined in eq. (3J.) without 



any cut on the dilepton invariant mass and for three different cuts as defined in eq. ( 35 ) 



asymmetry were worked out in both the HQET and FM model approaches. We find that the 
hadron energy spectrum is stable against variation of the FM model parameters. However, 
the hadronic invariant mass is sensitive to the input parameters. This dependence was already 
studied for the SD-contribution in [|],||] both in the context of the FM model and HQET. 

• We have quantitatively studied the uncertainties related to the implementation of the resonant 

rr 

and non-resonant parts in the coefficient Cg (s) in B —* X s £ + £~ . The numerical differences 
between the approach followed in [ 10 1 and the alternative ones, discussed in jyj] and j^], are 
found to be small in the dilepton invariant mass spectrum and negligible in the hadron energy 
and invariant mass spectra and spectral moments. In contrast, theoretical spectra are found to 
be more sensitive to the parameters Ai and A. 

• We have studied the hadron spectra by imposing the experimental cuts designed to suppress 
the resonant cc contributions, as well as the dominant BB background leading to the final 
state BB — ► X s £ + £~ (+ missing energy). In particular, the survival probability of the B — ► 
X s ^ + £~signal resulting from imposing a cut on the hadronic invariant mass Sh < 3.24 GeV 2 , 
as used in the CLEO analysis, is estimated and its parametric dependence studied. We have 
shown that the cuts such as the ones used in |7j effectively suppress the resonant contribution. 
Thus, the cut spectra essentially test the physics of the short-distance (and non-resonant cc) 
contribution, which can be systematically studied in perturbation theory and HQET. 

We hope that this work which provides a detailed theoretical profile of the hadron spectra in 
the decay B — > X s £ + £~w\W. be helpful in experimental searches of the rare decay B — > X s £ + £+. 
The distributions presented here will allow direct comparison of data with SM and will be useful in 
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FM 


No s-cut 


No s-cut 


cut A 


cut B 


cut C 


parameters 






e+ 


e 




1* 


e+ 


e 


£+£- 


(Ai,A) 


(Sh) 




(Sh) 


(Sh) 


(Sh) 


(S 2 H ) 


(Sh) 


(Sh) 


(Sh) 


(S 2 H ) 


GeV 2 , GeV 


GeV 


GeV 4 


GeV 


GeV 4 


GeV 


GeV 4 


GeV 2 


GeV 4 


GeV 


GeV 4 


(-0.3,0.5) 


1.47 


2.87 


1.52 


3.05 


1.62 


3.37 


1.66 


3.48 


0.74 


0.69 


(-0.1,0.4) 


1.57 


2.98 


1.69 


3.37 


1.80 


3.71 


1.88 


3.99 


0.74 


0.63 


(-0.14,0.35) 


1.31 


2.34 


1.38 


2.55 


1.47 


2.83 


1.52 


2.97 


0.66 


0.54 


(-0.3,0.5) tot 


1.41 


2.61 


1.41 


2.62 


1.61 


3.32 


1.66 


3.47 


0.74 


0.68 


(-0.1,0.4)** 


1.35 


2.14 


1.36 


2.15 


1.77 


3.60 


1.87 


3.94 


0.74 


0.62 


(-0.14,0.35)^ 


1.17 


1.84 


1.18 


1.85 


1.45 


2.76 


1.51 


2.95 


0.66 


0.54 


A 


0.15% 


0.19% 


0.22% 


0.42% 


0.90% 


1.56% 


0.32% 


0.58% 


0.01% 


0.32% 



Table 6: Spectral moments {Sh) and (S H ) for B — ► X s £ + £ , £ = fi,e for different FM model parame- 
ters and a hadronic invariant mass cutS H < 3.24 GeV 2 are given in the second to fifth columns. The 
values in the sixth to eleventh columns have additional cuts on the dilepton invariant mass spectrum 
as defined in eq. (35). The Sn-moments with cuts are defined in eq. fisHj). Entries in the first three 
rows are calculated using the SD- contribution alone. The subscript tot = SD + LD denotes that both 
the short and the long distance contribution are included in these moments. The values of A given 
in the last row represent the maximum uncertainty on the spectral moments (in%) resulting from the 
three approaches to take into account the continuum/ cc-resonant contributions discussed in the text. 



estimating the effects of various experimental cuts on the hadronic and dilepton invariant masses and 
hadron energy, which will be invoked in experimental analyses. Finally, this work underscores the 
importance of systematically improving the theoretical precision of the hadron spectra and spectral 
moments in the SD-contribution in the decay B — > V s ^ + ^ _ beyond what has been already done in 
as the theoretical uncertainties from the LD-contributions can be brought under control by judicious 
experimental cuts. 
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Figure 7: Hadron energy spectrum in B — > X s £ + £ in the Fermi motion model with the cuts on the 



dilepton mass defined in eq. (35); (a),(c),(e) without and (b),(d),(f) with the cc-resonance contribu- 
tion corresponding to cut A,B,C, respectively. The solid, dotted and dashed curves correspond to the 
parameters (Ai,A) = (-0.3, 0.5), (-0.1, 0.4), (-0.15, 0.35) in (GeV 2 , GeV), respectively. 



23 



(a) 



(b) 




(e) 



(f) 



4 6 

S M [GeV 2 ] 



S M [GeV 2 ] 



Figure 8: Hadronic invariant mass spectrum in B — > X s £ + £~in the Fermi motion model with the 
cuts on the dilepton mass defined in eq. $3?\); (a),(c),(e) without and (b),(d),(f) with the cc-resonance 
contribution corresponding to cut A,B,C, respectively. The solid, dotted and dashed curves correspond 
to the parameters (Ai,A) = (-0.3, 0.5), (-0.1, 0.4), (-0.15, 0.35) in (GeV 2 , GeV), respectively. 
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